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One of the most important requirements for a detector at the ILC is good jet energy resolution. 
It is widely believed that the particle flow approach to calorimetry is the key to achieving the 
ILC goal of a di-jet invariant mass resolution Om/m < Fz/mz. This paper describes the current 
performance of the PandoraPFA particle flow algorithm. For simulated light quark jets in the 
Tesla TDR detector, the jet energy resolution achieved is better than as /E ~ 3.4 % for jet energies 
in the range 45 — 250 GeV. This represents the flrst demonstration that Particle Flow Calorimetry 
can reach the ILC jet energy resolution goals. 
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I. INTRODUCTION 

^ ■ Many of the interesting physics processes at the ILC will be characterised by multi-jet final states, often 

I ' accompanied by charged leptons and/or missing transverse energy associated with neutrinos or the lightest 
CZ2 , super-symmetric particles. The reconstruction of the invariant masses of two or more jets will provide a powerful 
^ ' tool for event reconstruction and identification. Unlike at LEP, where kinematic fitting[i[ enabled precise 
jet-jet invariant mass reconstruction almost independent of the jet energy resolution, at the ILC this mass 
Q ■ reconstruction will rely on the detector having excellent jet energy resolution. The ILC goal is to achieve a mass 
I resolution for W q'q and Z — > qq decays which is comparable to their natural widths, i.e. am/m — 2.7% « 
J>-j' Tw/rnw ~ T^/m^. For a traditional calorimetric approach, a jet energy resolution of <te/E = ctj \/ i?(GeV) 
. leads to a di-jet mass resolution of roughly Um/m = a/ Ejj{GeV), where Ejj is the energy of the di-jet 
system. At the ILC typical di-jet energies will be in the range 150 — 350 GeV, suggesting the goal of as/E ~ 
0.3/ y^E{GeV). This is more than a factor two better than the best jet energy resolution achieved at LEP, 
<te/E = 0.6(1 -t- |cos6i|)/v/£^(GeV) |2|]. Meeting the jet energy resolution goal is a major factor in the overall 
design of a detector for the ILC. 
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II. THE PARTICLE FLOW APPROACH TO CALORIMETRY 



0^ 

' It is widely believed that the most promising strategy for achieving the ILC jet energy goal is the particle 
C , flow analysis (PEA) approach to calorimetry. In contrast to a purely calorimetric measurement, PEA requires 
' the reconstruction of the four-vectors of all visible particles in an event. The reconstructed jet energy is the 
^ , sum of the energies of the individual particles. The momenta of charged particles are measured in the tracking 
" '~j ' detectors, while the energy measurements for photons and neutral hadrons are obtained from the calorimeters. 
r> I The crucial step in PEA is to assign the correct calorimeter hits to reconstructed particles, requiring efficient 
^ . separation of nearby showers. 
■ ■ ■ ' Measurements of jet fragmentation at LEP have provided detailed information on the particle composition 

of jets (e.g. 0, Q)- On average, after the decay of short-lived particles, roughly 62% of the energy of jets is 
carried by charged particles (mainly hadrons), around 27% by photons, about 10% by long-lived neutral hadrons 
{e.g. u/Kl), and around 1.5% by neutrinos. Assuming calorimeter resolutions of cte/E — 0.15/-\/i?(GeV) for 
photons and ge/E — 0.55\/i?(GeV) for hadrons, a jet energy resolution of 0.19/-\/i?(GeV) is obtained with 
the contributions from tracks, photons and neutral hadrons shown in Tab. |T1 In practise it is not possible to 
reach this level of performance for two main reasons. Firstly, particles travelling at small angles to the beam 
axis will not be detected. Secondly, and more importantly, it is not possible to perfectly associate all energy 
deposits with the correct particles. For example, if a photon is not resolved from a charged hadron shower, 
the photon energy is not counted. Similarly, if part of charged hadron shower is identified as a separate cluster 
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the energy is effectively double-counted. This confusion degrades particle flow performance. Because confusion, 
rather than calorimetric performance, determines the overall performance, the jet energy resolution achieved 
will not, in general, be of the form as/E = a/ y^E{GeV). 

The crucial aspect of particle flow is the ability to correctly assign calorimeter energy deposits to the correct 
reconstructed particles. This places stringent requirements on the granularity of electromagnetic and hadron 
calorimeters. Consequently, particle flow performance is one of the main factors driving the overall ILC detector 
design. It should be noted that the jet energy resolution obtained for a particular detector concept is the 
combination of the intrinsic detector performance and the performance of the PFA software. 



Component 



Detector Energy Fraction Energy Res. Jet Energy Res. 



Charged Particles (X='=) 


Tracker 


~ 0.6 Bjet 


10-* 




< 3.6 X 10-^ Ef^^ 


Photons (7) 


ECAL 


~ O.SSjet 


0.15 




0.08 ^/E~t 


Neutral Hadrons (ft") 


HCAL 


~ O.lSjet 


0.55 




0.17 



TABLE L Contributions from the different particle components to the jet-energy resolution (all energies in GeV). The 
table lists the approximate fractions of charged particles, photons and neutral hadrons in a jet and the assumed single 
particle energy resolution. 



III. THE PANDORAPFA PARTICLE FLOW ALGORITHM 

PandoraPFA is a C-|--f implementation of a PFA algorithm running in the Marlin[^, framework. It was 
designed to be sufficiently generic for ILC detector optimisation studies and was developed and optimised using 
events generated with the Mokka?| program, which provides a GEANT4'8l simulation of the Tesla TDR 9| 
detector concept. The PandoraPFA algorithm performs both calorimeter clustering and particle flow in eight 
main stages: 

i) Tracking: for the studies presented in this paper, the track pattern recognition is performed using Monte 
Carlo information^. The track parameters are extracted using a helical fit. The projections of tracks onto 
the front face of the electromagnetic calorimeter are calculated using helical fits (with no accounting for energy 
loss along the track). Neutral particle decays resulting in two charged particle tracks (F°s) are identified by 
searching from pairs of non- vertex tracks which are consistent with coming from a single point in the central 
tracking chamber. Kinked tracks from charged particle decays to a single charged particle and a number of 
neutrals are also identified, as are interactions in the tracking volume (prongs). 

ii) Calorimeter Hit Selection and Ordering: isolated hits, defined on the basis of proximity to other hits, 
are removed from the initial clustering stage. The remaining hits are ordered into pseudo-layers which follow 
the detector geometry so that particles propagating outward from the interaction region will cross successive 
pseudo-layers. The assignment of hits to pseudo-layers removes the dependence of the algorithm on the explicit 
detector geometry whilst following the actual geometry as closely as possible. Within each pseudo-layer hits 
are ordered by decreasing energy. 

iii) Clustering: the main clustering algorithm is a cone-based forward projective method working from inner- 
most to outermost pseudo-layer. In this manner hits are added to clusters or are used to seed new clusters. 
Throughout the clustering algorithm clusters are assigned a direction (or directions) in which they are growing. 
The algorithm starts by seeding clusters using the projections of reconstructed tracks onto the front face of the 
calorimeter. The initial direction of a track-seeded cluster is obtained from the track direction. The hits in 
each subsequent pseudo-layer are then looped over. Each hit, i, is compared to each clustered hit, j, in the 
previous layer. The vector displacement, ry, is used to calculate the parallel and perpendicular displacement 
of the hit with respect to the unit vector(s) u describing the cluster propagation direction(s), = rij.u and 
d± = |rij X u|. Associations are made using a cone-cut, d± < d\\ tana -|- (3Dpad, where a is the cone half-angle, 
Dpad is the size of a sensor pixel in the layer being considered, and /3 is the number of pixels added to the 
cone radius. Different values of a and /3 are used for the ECAL and HCAL with the default values set to 
{tanaE = 0.3, /9e = 1-5}, and {tan an = 0.5, /3h = 2.5} respectively. Associations may be made with hits in the 
previous 3 layers. If no association is made, the hit is used to seed a new cluster. This procedure is repeated 
sequentially for the hits in each pseudo-layer (working outward from ECAL front-face). 

iv) Topological Cluster Merging: by design the initial clustering errs on the side of splitting up true 
clusters rather than clustering energy deposits from more than one particle. The next stage of the algorithm is 
to merge clusters from tracks and hadronic showers which show clear topological signatures of being associated. 
A number of track-like and shower-like topologies are searched for including looping minimum ionising tracks. 
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back-scattered tracks and showers associated with a hadronic interaction. Before clusters are merged, a simple 
cut-based photon identification procedure is applied. The cluster merging algorithms are only applied to clusters 
which have not been identified as photons. 

v) Statistical Re-clustering: The previous four stages of the algorithm were found to perform well for jets 
with energy less than 50 GeV. However, at higher energies the performance degrades rapidly due to the 
increasing overlap between hadronic showers from different particles. To address this, temporary associations 
of tracks with reconstructed calorimeter clusters are made. If the track momentum is incompatible with the 
energy of the associated cluster re-clustering is performed. If -Ecal — -Etrack > S.Sctb, where aE is the energy 
resolution of the cluster, the clustering algorithm, described in Hi) and iv) above, is reapplied to the hits in 
that cluster. This is repeated, using successively smaller values of the as and /3s in the clustering finding 
algorithm (stage Hi)) until the cluster splits to give an acceptable track-cluster energy match. Similarly, if 
£'TRACK — EcAh > 3.5tT£; the algorithm attempts to merge additional clusters with the cluster associated with 
the track. In doing so high energy clusters may be split as above. 

vi) Photon Recovery and Identification: A more sophisticated photon identification algorithm is then 
applied to the clusters. The longitudinal profile of the energy deposition, AE, as a function of number of 
radiation lengths from the shower start, t, is compared to that expected for an electromagnetic shower: 

AEk. Eo ^ ' ' , , At where a = 1.25 + -In—, 

r(a) 2 Ec 

Eq is the shower energy and E^ is the critical energy which is tuned give the appropriate average shower profile 
in the ECAL. The resulting level of agreement is used to improve the tagging of photons and to recover primary 
photons merged with hadronic showers. 

vii) Fragment Removal: At this stage there is still a significant number of "neutral clusters" (not identified 
as photons) which are fragments of charged particle hadronic showers. An attempt is made to identify these 
clusters and merge them with the appropriate parent cluster. All non-photon neutral clusters, i, are compared 
to all clusters with associated tracks, j. For each combination a quantity, e^j , is defined which encapsulates the 
evidence that cluster i is a fragment from cluster j. The requirement, Rij, for the clusters to be merged, i.e. the 
cut on eij, depends on the location of the neutral cluster and the change in the for the track— cluster energy 
consistency that would occur if the clusters were merged, Ax^ = {Et^ack — Ej)"^ /a'^ — {E^^iACK — Ei — Ej)"^ /a'^. 
If eij > Rij the clusters are merged. This ad hoc procedure gives extra weight to potential cluster matches 
where the consistency of the track momentum and associated cluster energy improves as a result of the match. 

viii) Formation of Particle Flow Objects: The final stage of the algorithm is to create Particle Flow 
Objects (PFOs) from the results of the clustering. Tracks are matched to clusters on the basis of the distance 
closest approach of the track projection into the first 10 layers of the calorimeter. If a hit is found within 
50 mm of the track extrapolation an association is made. If an identified kink is consistent with being from a 

fx^v or fi^v decay the parent track is used to form the PFO. The reconstructed PFOs are written 

out in LcioQ format. 

IV. CURRENT PERFORMANCE 

Fig. [T^) shows an example of a PandoraPFA reconstruction of a 100 GeV jet from a Z — > uu decay at 
\/s — 200 GeV. The ability to track particles in the high granularity Tesla TDR calorimeter can be seen clearly. 
Fig. lb) shows the total PFA reconstructed energy for Z uds events with |cos6'qq| < 0.7, where 6'qq is the polar 
angle of the generated qq system. These events were generated at ^/s = 91.2 GeV using the Tesla TDR detector 
model with a HCAL consisting of 63 layers and in total 6.9 interaction lengths. The root-mean-square deviation 
from the mean (rms) of the distribution is 2.8 GeV. However, quoting the rms as a measure of the performance 
over-emphasises the importance of the tails. It is conventional to quote the performance in terms of of rmsgo, 
which is defined as the rms in the smallest range of reconstructed energy which con tains 90 % of the events. 
For the data shown in Fig. (Tja) the resolution achieved is rmsgo/i? = 0.23/ ^^^(GeV), equivalent to a single jet 
energy resolution of 3.3%. The majority of interesting ILC physics will consist of final states with at least six 
fermions, setting a "typical" energy scale for ILC jets as approximately 85 GeV and 170 GeV at ^/s = 500 GeV 
and y/s =lTeV respectively. Fig. [2] shows the jet energy resolution for Z ^uds events plotted against |cos6'qq| 
for four different values of -y/s. The current performance is summarised in Tab. [Tll The observed jet energy 
resolution in simulated events is not described by the expression cje/E = a/ E{Ge\f). This is not surprising, 
as the particle density increases it becomes harder to correctly associate the calorimetric energy deposits to the 
particles and the confusion term increases. The table also shows a measure of the single jet energy resolution, 
obtained by dividing rmsgo by ^/2. For the jet energies considered (45 — 250 GeV) the fractional energy resolution 
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FIG. 1; a) PandoraPFA reconstruction of a 100 GeV jet in the Mokka simulation of the Tesla TDR detector, b) The 
total reconstructed energy from reconstructed PFOs in Z ^ uds events for initial quark directions within the polar angle 
acceptance jcos^qqj < 0.7. The solid line shows a fit to two Gaussians with a common mean; the broader Gaussian is 
constrained to contain 25 % of the events. The narrow Gaussian has a width of 2.2 GeV. 



is significantly better than the ILC requirement of 3.8% obtained from the consideration of gauge boson di-jet 
mass resolution. It should be noted that in a real physics analysis the performance is likely to be degraded by 
jet finding, jet-pairing and the presence of missing energy from semi-leptonic heavy quark decays. Nevertheless 
the results presented in this paper already provide a strong indication that Particle Flow Calorimetry will be 
able to deliver the ILC jet energy goals and it is expected that the performance of PandoraPFA will improve 
with future refinements to the algorithm. 



V. CONCLUSIONS 



Particle flow calorimetry is widely believed to be the key to reaching the ILC jet energy resolution goal of 
a di-boson mass resolution of CTm/m < 2.7%. Consequently, the design and optimisation of detectors for the 
ILC depends both on hardware and on sophisticated software reconstruction. Based on the PandoraPFA 
reconstruction of simulated events in Tesla TDR detector concept, it has now been demonstrated that particle 
flow calorimetry can meet this goal at the ILC. This was not true at the time of LCWS06 and, thus, represents 
a significant step forward in the design and future optimisation of the ILC detector (s). 



Jet Energy 


rmsgo 


rmS9o/A/£^jj(GcV) 


rmsgo / \/2Ej 


45 GeV 


2.2 GeV 


23% 


3.3% 


100 GeV 


4.1 GeV 


29% 


2.9% 


180 GeV 


7.4 GeV 


39% 


2.9% 


250 GeV 


12.0 GeV 


54% 


3.4% 



TABLE II: Jet energy resolution for Z — >uds events with |cos6'qqj < 0.7, expressed as, rmsgo for the di-jet energy 
distribution, the effective constant a in rmsgo/iS = a{Ejj)/ \J Ejj{GeY), and the fractional jet energy resolution for a 
single jet. 
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FIG. 2: The jet energy resolution, defined as the a in as/E — a\/E{GeV), plotted versus cosdqq for four different 
values of ^/s. 
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